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Inflation is the leading paradigm for explaining the origin of primordial density perturbations and 
the observed temperature fluctuations of the cosmic microwave background. However many open 
questions remain, in particular whether one or more scalar fields were present during inflation and 
how they contributed to the primordial density perturbation. We propose a new observational test 
of whether multiple fields, or only one (not necessarily the inflaton) generated the perturbations. We 
show that our test, relating the bispectrum and trispectrum, is protected against loop corrections 
at all orders, unlike previous relations. 



Introduction: The statistical distribution of the primordial density field provides a unique opportunity to test our 
understanding of the origin of the observed Universe. Inflation is the leading paradigm, but many questions about its 
details remain and it is so flexible that it is unclear if it can ever be ruled out. However it is at least possible to test, 

■ in a model independent way, whether a single source present during inflation was responsible for generating density 
perturbations, or whether multiple sources are required, by exploiting deviations from a Gaussian distribution. Non- 
Gaussianity (nonG) parameters are related by distinctive consistency relations, whose structure generally depends on 
whether perturbations are produced by one or more scalar fields. These parameters could be measured by forthcoming 
Planck satellite data and we may be able to answer fundamental questions about the number of degrees of freedom 

^) • contributing to the primordial density perturbation. 

, In single-source scenarios, density perturbations are generated by quantum fluctuations in a single scalar field, 

■ that does not necessarily correspond to an inflaton field driving the Hubble expansion. Single-field inflation is one 
particular case of this set-up, in which a single inflaton field both drives inflation and generates primordial cosmological 

■ perturbations. More generally, we include also set-ups such as curvaton or modulated reheating scenarios, in the limit 
J> 1 that only one field generates the density perturbation [1]. In these models, large nonG of local shape can be induced by 

a single light scalar field, whose dynamics becomes important only after inflation ends. In multiple-source scenarios, 
instead, the density perturbation is produced by more than one scalar field. This category includes models of multi- 
field inflation and multiple-field versions of other scenarios; all of which may also produce large local nonG. 

One famous consistency relation associates the squeezed limit of 3-point function with the scale dependence of the 
2-point function: /nl = — 5 / 12 ( n C ~ -0- This relation is valid only for single-field (clock) inflation [2], and is violated 
in more general single-source, or multiple-source scenarios. Another consistency relation, on which we will focus, 
connects the collapsed limit of 4-point function to the squeezed limit of the 3-point function: 



TNL = ( 6 W5) 2 . (1) 



This equality is satisfied at tree-level in general single-source scenarios, while it is generally violated in multiple-source 
set-ups, leading to the more general Suyama-Yamaguchi inequality tnl > ( 6 ^ NL / 5 ) [3-6]. However, as recently pointed 
out [7, 8], the equality (1) can be broken at an observable level, even in single-source scenarios, due to loop corrections. 

Inflationary observables associated with any given n-point (n-pt) function receive loop contributions, in terms of 
integrals over internal soft momenta that induce logarithmic corrections proportional to parameters related with higher 
n-pt functions. These contributions are particularly clearly seen in a diagrammatic representation of n-pt functions 
in terms of Feynman-type diagrams [9]. Loop corrections to correlation functions contribute to /nl and tnl, and 
in single-source scenarios, these corrections can combine in such a way to apparently break the equality (1). As we 
will discuss, the violation of (1) can be interpreted as due to the fact that the equality written in this way does not 
include contributions of soft momentum lines, connecting the different n-pt functions entering the definitions of the 
parameters tnl and /^l- These soft momentum lines are allowed by momentum conservation, and lead to radiative 
corrections comparable to the conventional loop corrections in the limits of collapsed and squeezed configurations. 
Accounting for both loop corrections and soft modes, we will build a new combination of bispectrum, trispectrum 
and power spectrum, that leads to an equality satisfied to all orders in radiative corrections in generic single-source 
scenarios. The equality reduces to equation (1) at tree-level, and is in general broken in multiple-field scenarios. Our 
result therefore generalizes the Suyama-Yamaguchi relation to all orders in radiative corrections. Moreover, in the 
second part of this work, we will also discuss how the soft modes we consider are the source of the inhomogeneity of 
nonG observables discussed in [7], leading to further observational implications for our findings. 

The role of soft momenta is reminiscent of what happens in QED, in which a careful inclusion of contributions of 
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both real and virtual soft photons is crucial for canceling IR divergences in physical processes [10] . The conceptual idea 
we develop here is similar to what happens in that context, although the technical implementation will be different. 

Radiative corrections to n-pt functions. We focus on a local Ansatz for the primordial curvature perturbation 

[11] 

C(x) = Cg(x) + — — CgW + ^-C G ( X ) + "j^— £g( x ) + (Cm) . ( 2 ) 

where Cg is a Gaussian random fluctuation, with vanishing ensemble average ((c) = 0. The local Ansatz assumes 
the parameters /nl, 3nl and /inl are constants. We use a bar to emphasize that these parameters, as defined by (2), 
characterize tree-level amplitudes of 3-, 4-, and 5-pt functions respectively. We focus on single-source scenarios, in 
which the curvature perturbation is generated by a single scalar a and with observably large non-Gaussianity. It is 
then possible to rewrite the local Ansatz (2) in terms of an expansion of a suitable classical function N(a), 

C(x) = N(a + 5a(x)) - (N(a + <5a(x))) = ]T ^ (fc(x)" - (5a(x) n )) (3) 

n— 1 

This expansion is similar to the <5A~ approach of [12] but for our purposes it is not necessary to specify N(a) further. 
o o denotes a homogeneous, time dependent background solution and the fluctuations 8a(x) are Gaussian with zero 
mean. Making the identification N'(<r )Sa(x) — Cg( x ), a comparison of the equations (2) and (3) yields the relations 

6/nl/ 5 = W'/at' 2 , 54g NL / 25 = JV'"/jV' 3 , etc [12]. 

The general definitions for the parameters /nl and tnl at tree-level and beyond are given in appropriate squeezed 
and collapsed limits by 

/nl = 12 a p (fcl) p(fc 2 ) ' (4) 

i ,. r(fci,fc 2 ,fc 3 ,fc4,fci2,fci 3 ) 

TNL = I fc !™o P(fc 12 )P(fci)P(fc 3 ) ' ( ' 

where (C kl Ck 2 ) - (2 7 r) 3 P(5(k 1 +k 2 ), (C kl Ck 2 Ck 3 ) = (2^) 3 P(5(k 1 +k 2 +k 3 ), (Ck, Ck 2 Ck 3 Ck 4 ) = (2^) 3 T(5(k 1 +k 2 +k 3 +k 4 ), 
and kij = |kj + kj|. Notice that at tree-level /nl defined by equation (4) reduces to the parameter /nl entering the 
local Ansatz (2). At tree level one has /nl = bN " /&n' 2 , tnl = n " 2 /n" l , from which equality (1) follows immediately. 
However, when 1-loop contributions are added to the tree-level results, one finds [15] 



/nl' P /m. ^ (2/nl - 3/NLffNL - 3%l) Pin ('^) , ((».) 
324 „_ 9 „ T y x ^, / fc ^ 



Tn° l op = tnl - ^5 ( 8 /nl ~ 12/nl.9nl - 9g NL - 12/inl/nl) Pin . ; , 



where 27r 2 P = k 3 P(k), and we neglect its weak scale dependence. Hence 

■Pin 



/ loop \ ^ 

loop 6/n°l° P 
NL - E 



25/ N 2 L \kmj 



(7) 



showing that the consistency relation (1) is violated already at one loop in single-source models, if the tree-level 
t/NL is non- vanishing. If the non-Gaussian parameters are large, this violation of the consistency relation can be 
observed by the Planck satellite. On the other hand, while being perfectly well-defined and consistent, the standard 
way of calculating loop corrections individually to the power spectrum, bispectrum and trispectrum, and then taking 
appropriate ratios to define the non-linearity parameters beyond tree-level, can miss important physical contributions. 
Indeed, these loop corrected quantities correctly characterize the ratios of individually measured bispectrum, or 
trispecturm, and power spectrum. However, when simultaneously measuring combinations of n-pt functions, such as 
the ratios in (4) and (5), one should allow for the inclusion of soft lines connecting distinct n-pt functions. Although 
momentum is of course conserved, no detector is sensitive enough to probe these soft lines: their contribution is 
physical and must be included. This observation suggests that, besides considering the relation (1), by including the 
effects of soft modes it is possible to build a new observable combination of n-pt functions that leads to an equality 
protected against radiative corrections. See also [13] for the slow-roll suppressed effect of soft modes on the power 
spectrum in single clock inflation. 

Diagrammatic approach to loop corrections. It is illuminating to discuss the role of soft modes diagrammatically, 
first in a simple example, and then applying our observations to equality (1). We implement the diagrammatic 
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approach of [9]. We use solid dots to mark external momenta k^, with the number of attached propagators to each 
vertex (corresponding to P) giving the number of derivatives of the function N, defined by equation (3), associated 
to the vertex. There are no internal vertices since we assume that 5a is Gaussian. The numerical factors are the total 
for each diagram relative to the tree- level term (which may have some possible permutations), and are given by the 
numerical factor of the given diagram ( 1 /2 if a dressed vertex, otherwise unity at one loop or tree- level), times the 
number of distinct ways in which the loop can be drawn onto the tree-level diagram. Up to one loop (the first diagram 
is always the unique tree-level diagram) the power spectrum, bispectrum, and the tnl diagram for the trispectrum 
are respectively given by 

— + 1^ (8) 

•♦•+ | + 2 + (9) 

The above method correctly calculates the loop corrections to each individual n-pt correlator. However when consider- 
ing a combination of correlators with soft internal momenta, as is for example the case with P 2 in the denominator of 
/nLj onc should take into account the contributions of soft momenta connecting different elements of the combination. 

Let us illustrate the role of soft modes, by considering radiative corrections to the square of the power spectrum as an 
example. We denote with (. . . ) ra( j the sum of tree-level and radiative contributions to a given quantity. We estimate, 
to first order in radiative corrections, the difference between (-P) 2 ad , associated with the observable (£(k)£(k')} 2 , and 
(~P 2 )rad, denoting another observable, (C 2 (k)£ 2 (k')). Diagrammatically: 



(P(fc)) r 2 ad - (P(fc) 2 ) rad 



O 2 ^^^^ A T"" 

rCH + A ) -(_+- 1^+2^+41+ 



2 ) V" 2 

= -4(6/5/ NL ) 2 P(fc) 2 Pln(Vfc^) . (11) 



The final diagram represents a 4-pt function with a soft internal line (the thick-red line) connecting two N" vertices 
(which can be done in four possible ways). This contribution can not be observationally distinguished from the 
product of two disconnected 2-pt functions and it must therefore be included in (P 2 ) ra d- So, accounting for loop 
corrections only, would not give the correct result for the observable (11), which demonstrates the importance of the 
soft modes. 

Contributions to radiative corrections associated with soft lines connecting different diagrams play a crucial role for 
characterizing equality (1) beyond tree- level. In order to analyze it diagrammatically, it is convenient to re-express 
it in terms of square of the bispectrum, and the product of the trispectrum with the power spectrum. At tree-level, 
(1) reads B 2 kl ^ = TP(k 5 ) k 12 ^o,k 5 ^o- Including loop corrections to P, B and T, see Eqs. (8-10), in appropriate 
squeezed and collapsed limits, the relation is 

(P) r 2 ad - (P) rad (T) rad = -4fl»— • = -4 (54/25 5nl) 2 P(k) 2 P^ PH k /^n) ■ (12) 

So, the equality is broken by a term proportional to tree-level ()nl, as discussed in the previous section. On the other 
hand, a straightforward calculation shows that the relation 

/o2\rad / pT->\rad 



+ ( 6 /s) 2 ( 3 W/nl + 9 5nl/ 4 /5 l + 15g NL + 6/^ L )pln(fe/fcz«)] , (13) 



instead leads to an equality that is preserved by radiative corrections. When evaluating the previous quantities, it is 
crucial to include diagrams describing contributions of soft modes connecting different elements of each combination 
(drawn in red on the diagrams below), respectively the 3 diagrams for (-B 2 ) rad and 2 diagrams for (PT) ra d. 



(14) 



When including these diagrams with their associated numerical factors (following the rules developed in [9]), one 
fulfills the equality (13), (P 2 ) ra d = (PP)rad, that generalizes the equality (1) to first order in radiative corrections. 
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So the key observation is that, in order to define an equality that remains valid when radiative corrections are included, 
one should consider radiative contributions to the entire combination of (-B 2 ) ra d and {PT) Ta d, that include both loop 
corrections and contributions from soft modes connecting different diagrams. The new combination (13) can be 
considered as a new inflationary observable: when going beyond tree level in a loop expansion, it allows to probe 
nonG parameters in a different way with respect to relation (1). 

An alternative approach to radiative corrections. Let us now reconsider the problem from a different 
point of view, which allows a straightforward generalization of our results to all orders in radiative corrections, and 
emphasizes the connection to inhomogeneities of non-Gaussian observables [7]. First, some preliminary material. 
Consider a wavenumber Q < aH which defines a length scale smaller than the observed universe. The fluctuations in 
(3) can be divided into long wavelength (LW) and short wavelength (SW) components with respect to this scale as 

Jq>Q ( 27r ) Jq<Q ( 27r ) 

ee 8a s {x) + Sa L {x.) . (15) 

For a Gaussian field, the LW fluctuations are uncorrelated with the SW modes, (5ox(k)<5er s (k')) = 0, and they both 
have vanishing ensemble averages, {8a l) = {8a s ) — 0. Up to cosmic variance, the ensemble averages correspond to 
spatial averages over the full observable sky. 

For measurements probing wavenumbers k > Q, or equivalently regions of size smaller than 1 /q, the LW modes 
8a l act as an approximately homogeneous background for 8a s . Indeed, the average of 8a l computed over a spherical 
patch of volume Vx , with the origin located at a fiducial point x , is given by [7] 

{Sa L {x)) Xo ~ 8a L {x ) ■ (16) 

The average (<5<7j,(x)) Xo consequently depends on the location xo of the patch. In general, the contribution of long- 
wavelength modes differs from patch to patch, generating variations in observables evaluated in different subhorizon 
patches across the sky. This has a cumulative effect, leading to a log-enhanced variance of the long-wavelength 
fluctuations over the entire observable universe 

= / 10? PM ^° '"(£;) • ,17 » 

ki R <q<Q 

where V a is the dimensionless spectrum of the fluctuations 8a. 

These observations provide a convenient way to compute radiative corrections to inflationary observables. Using 
equations (3), (15) and (16), the curvature perturbation as measured within a patch of volume V Xo can be written as 

x ° ' n\ ' n! 

n— 1 n— 1 

Converting this expression to Fourier space on scales k > Q is straightforward, as <5<tl(xo) acts as a constant under 
this operation. On the other hand, the LW fluctuations are operators with non-vanishing ensemble, or full sky, 2- 
pt functions (17). This has interesting consequences when considering full sky expectation values of tree- level n-pt 
functions evaluated in small patches ( (£(ki)...£(k„)) Xo } . We can expand the argument of (...) in terms of LW modes 
8a l , and evaluate the ensemble averages: at this point the connection with loops in inflation becomes apparent. 
Indeed, this operation is equivalent to computing radiative corrections to the corresponding full sky n-pt functions 
in a leading-log approximation. Ensemble averages of powers of 8a l lead to log-enhanced contributions, controlled 
by formula (17). By identifying the reference scale Q with the scale k at which the measurement of full sky n-pt 
functions are performed, one finds that the LW mode contributions exactly reproduce the radiative corrections in the 
leading-log approximation. 

More specifically: both radiative corrections and LW contributions arise from the expansion of the curvature 
perturbation (18) around the same classical background solution. The vertices appearing in the loop diagrams, and 
diagrams containing soft modes are by construction identical to the ones associated with LW corrections. Integrals 
over the loop momenta can be estimated using the so called pole approximation [14], which amounts to integrating 
each of the loop propagators from the IR cutoff fciR. up to the mode k at which the corrections are computed. This 
produces n logarithms of the form ln( fc /fe m ) where n is the number of lines where the loop momentum appears. 
Exactly the same n logarithms are reproduced by expanding the same correlator in long-wavelength modes provided 
that we set Q — > k. The same applies for integrals over soft modes. Indeed, since the radiative corrections consist 
of long-wavelength modes, and since we are performing a semiclassical leading logarithm expansion of the curvature 
perturbation, the exact correspondence is precisely what one expects to find. 
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This approach allows us to easily reproduce and extend our discussion of radiative corrections to the equality (1). 
From the expansion (18) we find the observables and fj^, measuring respectively the squeezed and collapsed 
limits of tree-level 3- and 4-pt functions within a small patch V Xo , are given by 



_B = 5 N"(a +6(T L (xo)) 

12 feV^o P(k 1 )P(k 2 ) 6 N' 2 (<j + 6a L {x )) 



Xo 1 f = N" 2 (a + Sa L (^)) 

NL 4 fc 12 ^o P(A5i2)P(fci)P(fc 3 ) ^ ,4 (ao + 5a L (x )) ' 

Taking ensemble averages of (/j^) 2 and we then obtain the equality 

<^> = ( 6 /5) 2 <(/£l) 2 > . (20) 
that can be expressed also as (-B XOiSq ) = (T XOiCO n Pxo.soft), corresponding to 

lim {Bl {k u k 2 ,k 3 )) = lim {T XQ {k u k 2 ,k 3 ,k 4 ,k 12 ,k 13 )P Xo (k 5 )) . (21) 

fei^O ki^-k 2 

&2— *-&3 fel— >fe;? 

fc 5 ->0 

As discussed above, this relation between tree-level quantities in a small patch V Xo , equates to all orders in radiative 

— 1/3 

corrections, and to leading logarithm accuracy, the corresponding full-sky quantities evaluated at the scale k ~ Vx . 
Expanding (20) to second order in the LW modes 8<jl exactly reproduces the first order radiatively corrected equality 
between (-B 2 ) ra d = (TP) la d, derived above using diagrammatic methods (13). In squaring one obtains quadratic 
contributions in So x both from the linear and quadratic terms in equation (19). The former correspond to soft 
corrections, and are the origin of the inhomogeneities of nonG discussed in [7], and the latter to loop corrections. 
Both of them have to be included in order for the equality (13) (or (20)) to be satisfied. This is why the the loop 
corrected quantities /nl P an d t nl P (6)> missing these soft contributions, in general fail to satisfy the equality (1). 
Focussing at one-loop order, it is instructive to rewrite (21) as 

(Bxo, sq ) 2 = ( T X0;Coll )< P xo , soft ) + yj (a 2 T a%) - 0% (22) 
= (^ co n)(P_ ft )(l-(^)VlnA). 

Here a 2 = (x 2 ) — (x) 2 denotes the non-zero variance generated by the long- wavelength modes, giving a statistical 
variation to quantities measured on subhorizon patches. Writing the formula in this form explicitly shows that taking 
averages over full sky of single quantities, and then combining them together, does not lead to a simple form of the 
equality. Additional pieces proportional to §nl (and /inl at two loops) would seem to lead to a violation of the tree- 
level result (1) when loops are included. Diagrammatically, as we have seen in the previous section, these corrections 
can be traced back to soft internal modes connecting tree-level correlators. This conclusion is avoided by properly 
defining the quantities that, once averaged over the full sky, allow us to write an equality in a form that automatically 
handles radiative corrections at all orders, as given by Eq. (21). 

Conclusions and Outlook: In this work we have investigated new contributions to non-Gaussian inflationary 
observables in squeezed or collapsed configurations, associated with soft momentum lines connecting different n-pt 
functions. Our analysis is essential for investigating and understanding consistency relations among inflationary 
observables, that can be tested by the Planck satellite and provide model independent information about the number 
of degrees of freedom contributing to the primordial curvature fluctuations. We showed that the new contributions 
we discussed are essential for defining a combination of power spectrum, bispectrum and trispectrum, equation (21), 
corresponding to a new equality among observables preserved by radiative corrections in single-source inflationary 

scenarios. This provides the natural generalization of the tree-level equality tnl = ( 6 / 5 /nl) , which is broken by 
loop corrections. We discussed our results adopting a convenient diagrammatic representation of inflationary n-pt 
functions in terms of Feynman diagrams. We also made a connection between these results and inhomogeneities of 
non-Gaussian observables, clarifying the relation between loop corrections and inhomogeneous non-Gaussianity. In 
order to do this, we employed a particularly simple method based on splitting long from short wave-length modes with 
respect to a fiducial scale, and exploited the fact that long wavelength mode contributions to inflationary observables 
behave in the same way as radiative corrections. 

In summary, we have shown from various points of view that long wavelength, soft modes can provide physical con- 
tributions to inflationary observables. While in this paper we focussed on single-source scenarios, it is straightforward 
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to generalize the method and our results to a generic multiple-field case. The tree-level equality (1) gets replaced 

by the inequality tnl ^ ( 6 / 5 /nl) [3]. According to our previous discussion, this translates into a new inequality 
between full-sky observables: 

(B 2 Y^ ^(PT)T^, k ^, (23) 

which holds to all orders in radiative corrections, and to leading logarithm precision. This inequality can be unam- 
biguously used to discriminate between single and multiple-source scenarios for generating primordial perturbations 
at arbitrary orders in radiative corrections, even when loops are included. The Suyama-Yamaguchi inequality and our 
inequality (23) above are two observable relations probing different physics when loops arc included, and capable of 
testing in different ways models which generate the curvature perturbation. We will explore this and other interesting 
issues elsewhere. 
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